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Abstract 

In this paper we try to introduce the ladder operators associated with the pseudo- 
harmonic oscillator, after solving the corresponding Schrodinger equation by using the 
factorization method. The obtained generalized raising and lowering operators naturally 
lead us to the Dirac representation space of the system which is very easier to work 
with, in comparison to the functional Hilbert space. The SU{1, 1) dynamical symmetry 
group associated with the considered system is exactly established through investigating 
the fact that the deduced operators satisfy appropriate commutation relations. This re- 
sult enables us to construct two important and distinct classes of Barut-Girardello and 
Gilmore-Perelomov coherent states associated with the system. Finally, their identities 
as the most important task are exactly resolved and some of their nonclassical properties 
are illustrated, numerically. 



1 Introduction 

Molecular vibrational and rotational spectroscopies are usually known as the important parts 
of molecular physics. The harmonic oscillator could be considered as a useful ground to explain 
the molecular vibrations, but this model is very ideal and simple and so has some particu- 
lar restrictions. Strictly speaking, the molecular vibrations are really anharmonic, so if their 
Schrodinger equation can be exactly solved, they may provide more suitable models for the 
description of vibrating molecules. As some examples of anharmonic oscillator molecular po- 
tentials we may refer to the Morse [1], Poschl- Teller [21 13] and Mie-type [1] potentials which are 
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three typical model potentials which the molecular vibrations can be described satisfactorily 
by means of them. In this work, we want to study another anharmonic potential, the so- 
called pseudoharmonic oscillator, which is the sum of the harmonic oscillator and the inversely 
quadratic potential [S], i.e., 

V(.) = -,n.V + --, (1) 

where m, uj and a respectively represent the mass of the particle, the frequency and the strength 
of the external field. Sometimes this system has been called isotonic potential [HI [TJ El [9] . In fact, 
this potential was proposed by Goldman and Krivchenkov [TU]. Some methods that have been 
used to solve the wave equation for the pseudoharmonic potential can be summarized as follows: 
Nikiforov-Uvarov method |Tll [12], polynomial solution [13], exact quantization rule IH] and 
Laplace transform approach [15]. Moreover, obtaining the exact solutions of the Schrodinger 
equation for the pseudoharmonic oscillator using a new method (and in general enlarging the 
circle of solvable quantum systems) is one of the main problems in quantum physics [T6] . 

In the present study we firstly want to solve the latter equation by the factorization method 
which simply arrives us to the exact solution of the eigenvalues and eigenfunctions of the 
considered system. Very recently this method is used for the construction of a new set of 
generalized coherent state and then the geometric phase of the states is investigated [17] (see 
also [m [m |20l [21] and references therein). Our motivation hes in the fact that, by using 
this method we will be able to find out the appropriate generalized creation and annihilation 
operators for the chosen system, by which we are simply lead to the Dirac space with its 
instruments (ket-bra space) provides us with a very simpler tool than the usual functional 
space. Consequently, we can construct two new, important and distinct classes of generalized 
coherent states associated with such physical system namely Barut-Girardello and Gilmore- 
Perelomov coherent states. After this exploration, due to the importance of the nonclassicality 
of quantum states in various theoretical and experimental fields of physics in recent decades 
[22| [23] we are naturally lead to the investigation of the nonclassically signs of our introduced 
states for completing the work. Moreover, prior to these examinations we will exactly establish 
the resolution of the identity of both classes of states, as the most important property of any 
generalized coherent state. 

Notice that our present work can be considered in the line of the dynamical generalization of 
coherent states [2H[25l[26]. Also, it is worth mentioning that, the obtained Barut-Girardello and 
Gilmore-Perelomov coherent states also can be really categorized in the "algebraic" [27 [ [28 [ [29] 
and "group theoretical" [30] [31] [32] generalization of coherent states, respectively. From the 
physical point of view, we would like to mention that the strong relation of the laser light 
to the well-known standard bosonic coherent states corresponds to harmonic oscillator is well 
established from the sixth decades ego by the pioneering works of Glauber, Sudarshan and 
Klauder. Indeed, this was only after this discovery that, hard works along their generalizations 
as well as physical realization were began. This is why their various manners of generalizations 
are yet of hot importance [SSj [Ml [351 ESI [37] . Henceforth, apart from presenting a powerful 
technique which can be used for other known potentials, due to the direct relation of our 
introduced coherent states to a particular physical potential (which as we mentioned in the 
above lines, may be more suitable for the description of real physical phenomena) we hope that 
our work suitably concerns to physical realization, too. 

As is mentioned before, this potential has received a lot of attention in recent years. In 
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addition, the used approach in this paper is conceptually near to the well-known SUSY quantum 
mechanics techniques which is an interesting approach for many authors [381 EH HOl SD Hll 
H3l m]. Altogether, it may clearly be seen that our method to the solution of the eigenvalue 
problem and the from the final obtained results are essentially new. In this way, to the best 
of our knowledge, the approach to the treatment of this particular problem, the corresponding 
ladder operators, the associated coherent states of the types which are obtained, in addition to 
establishing the resolution of the identity for both classes of states and the related nonclassicality 
signs have not been appeared in the literature up to now. 

The organization of the paper is as follows. In the next Section, we are going to review the 
factorization method briefly. This method is used to solve the Schrodinger equation for pseu- 
doharmonic oscillator in one-dimensional space in Section 3. According to our obtained results, 
we construct the creation and annihilation operators, by which we study the dynamical group 
corresponding to this system in Section 4. Then, after producing the Barut-Girardello and 
Gilmore-Perelomov coherent states associated with the pseudoharmonic potential, the resolu- 
tion of the identity corresponding to both classes of introduced states is established in Section 
5. We demonstrate some of the nonclassical properties of the states like squeezing and sub- 
Poissonian statistics, numerically in Section 6. Finally, summary and conclusion are presented 
in Section 7. 



2 The factorization method: a brief review 

The factorization method is a useful technique that reduces the dynamical equation of a given 
system into a simple one which is easier to handle. This method was first introduced by 
Schrodinger |45l |46l |47] and then followed by Dirac [IH] to avoid the usage of cumbersome 
mathematical tools and it has also been successfully applied to resolve some problems for 
which there exists exact solutions. This method further developed by Infeld and Hull [19] and 
Green [50]. The characteristics and advantages of the factorization method can be summarized 
as follows. First, this method applies only to the discrete energy spectra. Secondly, we can 
achieve immediately the desired eigenvalues and the normalized eigenf unctions from the given 
Hamiltonian without any need to use the traditional methods in quantum mechanics. Thirdly, 
it is possible to avoid deriving the normalization constant, which is sometimes difficult to obtain. 
The fourth property comes from the fact that we may discover the hidden symmetry of the 
quantum system through the establishment of an appropriate Lie algebra, which usually can 
be realized by means of the ladder operators. This method is a very powerful tool for solving 
second-order differential equations and has attracted much attention of many authors. In this 
method, which is closely related to SUSY quantum mechanics, the second-order differential 
operator Hq (Hamiltonian) is written as the product of two first-order differential operators oq 
and Oq plus a constant Eq, 

Ho = ajao + Eq. (2) 

Then, rewriting the SUSY Hamiltonians H_ and respectively as Hq and Hi, and changing 
suitably the subscripts we have Hi = aoal + Eo. Also, by factorizing Hi one has Hi = a\ai + Ei, 
so firstly we have 

Hi = aoo-o + Eq = a\ai + Ei, (3) 
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with the SUSY partner of Hi as H2 = aia[ + Ei. Continuing the above procedure leads to the 
generation of a hierarchy of Hamiltonians [511 [52] given by 

Hn+i = a|^+ia„+i + En+i = a^al^ + En, n = 0, 1, 2, 3, ... . (4) 

The complete set of eigenfunctions can be generated by the successive action of these operators 
on the key eigenfunctions, which are the exact solutions of the first-order differential equation. 
Then, the recurrence relation may be used to achieve all eigenvalues. 
In fact, with the following theorem we can find eigenfunctions and eigenvalues [53] . 
Theorem: Assume that equations (I2l) to (jl]) with real constants En are true. Also suppose 
that there exists a null eigenfunction |^„) with zero eigenvalue for each a„, i.e., 

an\in) = 0. (5) 

Then, 

a) En is the nth eigenvalue of the Hamiltonian. 

b) The corresponding eigenfunction of the Hamiltonian H is, 

iV'n) = c„ajal...4_j^„), (6) 
where Cn is a normalization constant may be determined as: 

\Cn\ = [{En - En-i){En " ^„-2)...(^n - EoT'^^ . (7) 

Three important worthwhile points lie in the theorem: 

i) This method applies only to the discrete energy spectra. 

ii) In the selection of a„, one should choose it in such a way that En gets the maximum 
value. 

iii) The theorem doesn't give us any information about the form of a„. Moreover, evidently 
the form of depends explicitly on the particular potential V{x). Anyway, since alao is equal 
to H except for a constant, thus must contain a linear momentum term to be consistent 
with the kinetic energy part of the Hamiltonian. Accordingly, a„ is assumed to have the form 

an = -^={P + ifn{x)), (8) 
V2m 

where p is the momentum operator and fn{x) is a real (and differentiable) unknown function 
of X. This function would be such that equation (jl]) is satisfied. 



3 Solving the Schrodinger equation of pseudoharmonic 
oscillator: factorization method 

Now, we are going to search for the eigenenergies and eigenfunctions of a particular system, i.e. 
pseudoharmonic potential in one-dimensional space, by the factorization method. Let us start 
with the Schrodinger equation for this potential, i.e., 

(^^ + Vix)^^ix) = E^ix), (9) 
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where the operator p is defined as 4^. According to the form of pseudoharmonic oscillator 
potential introduced in ^ and keeping in mind ([8]) we can guess that fn{x) can take the 
following form 



X 



(10) 



where hn and c„ are constant values which should be chosen such that En gets the maximum 
value. Thus, a„ can be written as: 



1 



2m 



p + ii hnX + — 

X 



With the help of a„ and its conjugate a\^, one can easily check that the following relations hold: 



dnCLn 



2m 

1 
2m 



p' + + hn i2Cn + h) + ^{Cn- K) 



X^ 



p" + hlx^ + hn {2Cn - n) + ^ (C„ + h) 



X" 



With substituting n = in f lT^ and using ^ it is required that 



1 

2m - 



Co 



p^ + hlx^ + 6o (2co + ^) + ^ (co - n) 



x^ 



1 2 2 h (y 

1 — muj X H 7^. 

2m 2 2m 



p 



The above equation will be satisfied, if 



2, ,2 



6q = m~uj 



-muj. 



co(co - ^) = s{s + 



Co 



(12) 
(13) 



(14) 

(15) 

(16) 
(17) 



6o(2co + h) + 2mEo = 0, 
where the allowed values of ho in f[T^ can satisfy the condition (ii) of the stated theorem for 
which Eq gets the maximum value. Also, note that the change of variable s = — ^ + \J a + \ in 

f|T6|) has been used. With substituting ([T5|) and f|T6|) in f|T7|) . the maximum of the ground state 
energy will be achieved as: 

^ . ^ 3^ 

Eq = fUjJ \ s + 



For attaining the other eigenvalues we can use the recurrence relation (j4]), from which one has: 



2m 



Cn+l 



p + hn+iX + hn+1 (2Cn+l + ^) H 5- (c„+i - /l) 



2m 



p^ + 6^a;^ + hn (2c„, - n) + ^ (c^ + /i) 



x^ 



+ En+l 
+ En. 



The relation f lT^ will be satisfied, if 



bl^. = hi 



Jn — UQ 



-mu, 



(19) 



(20) 
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Cn+l{Cn+l -h) = Cn{Cn + K) > C„ = (s + ?2 + l)h, (21) 

6„+i(2c„+i + h) + 2mE„+i = 6„(2c„ - K) + 2mE„. (22) 
Finally, from (120|) . ( 12T|) and ( 122|) we can derive the required eigenvalues as 

E„ = (^2n + s + , n = 0,1,2,.... (23) 

With substituting 6„ and c„, respectively from fl2U]) and fl2T]) in flTT]) . a„(a;) will be explicitly 
obtained 

./ is + n + lW 

p + I \ —mux H — 



an{x) 



2m 



X 



(24) 



The equation (l23i) gives all eigenvalues, if there exists an appropriate answer for any arbitrary 
n in ([5]). Thus, we need to solve the following equation 

an{x)Ux) = (25) 

or equivalently 

x] ^n{x). (26) 



dx \ X ^ , 

The solution of the above differential equation can be simply obtained as follows 

^„(a;) = x^+"+ie-^^'. (27) 

By using (|6]) and after doing some tedious calculations, the normalized eigenfunctions of pseu- 
doharmonic oscillator can be obtained as {h = m = u = 1): 

Mx) = Nnx'^'e-'TLn^Hx^), (28) 

where Ln^^(a;^) is the associated Laguerre polynomials and the normalization constant in f l28l) 
is determined as 



/ 2r(n + l, , , 

Nn = \ L, \ 29 
L (n + s + ^ 



2' 



4 Introducing ladder operators associated with pseudo- 
harmonic oscillator 

It is remarkable that the operators a„ in (12^ and its conjugate do not play the role of raising 
and lowering operators for the eigenstates of the pseudoharmonic potential. But, keeping in 
mind these operators, we want to construct a rather new pair of operators which yield our 
mentioned requirement. For this purpose, we introduce a new operator as follows 

d s + n + 1 , , 

An = --r-x + , 30 

dx X 
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where its conjugate read as 

Ai = — -x + . (31) 

By applying the operator An on ipn{x) and using the following recurrence relations for the 
associated Laguerre polynomials |5i 



dx^^ ' ~\{n + l)L-+i(a;) -{n + a + l- x)L^{x), ^^^^ 



we have 



n 2 1 

Anipnix) = tpnix) + - \ n{n + s + -) ^Jn-i^x). (33) 

X X \ I 



Similarly, by acting A^^ on '^nix) one has: 



n + 1 2 / 3 

Ai^n{x) = ^nix) + -d{n + l){n + s + -) ^„+i(x). (34) 

Clearly, the above operators A^-, A}^ do not yet possess the role of ladder operators. Alto- 
gether, it is possible to find out the ladder operators which can regularly raise and lower any 
eigenf unctions of our considered system, i.e., 

_ _ 2/ d _ _ s 1 oc , , 

where we have introduced the number operator which satisfies the eigenvalue equation 

N\n) = n\n), ipn{x) = {x\n) (36) 

with {\n)}'^^Q as the eigenkets of the pseudoharmonic oscillator system. Notice that 7^ 
M+M_. It is easy to check that the ladder operators M± have the following properties 

M±\n) = m±\n ± 1) , n = 0,1,2,... (37) 

where 

m^ = ^in + -±-){n + s + -^). (38) 

We now study the dynamical group associated with the operators M+ and M„. Based on the 
obtained results in (1371) and (l38i) . we can calculate the following commutator: 

[M_,M+]\n) = 2mo\n), (39) 

where we have defined 

mo = n + I + -. (40) 

Thus, we can define the operator 

Mo = N + ^ + ^. (41) 
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Adding our above results, the operators M± and Mq justify the following algebraic relations 

[M_, M+] = 2Mo, [Mo, M±] = ±M±. (42) 

As is well known, according to f H2|) . the operators M-t and Mq correspond to an sm(1, 1) Lie 
algebra. In addition, the system Hamiltonian can be simply expressed as: 

H\n) = {2n + s + ^)\n) = [M_,M+]|n). (43) 

The form of the Hamiltonian in fHSl) satisfy the eigenvalue equation H\n) = En\n) with En de- 
termined in ( 123|) . It may be interesting to notice that while for the harmonic oscillator one has 
HHo\n)HO = {n + l/2)\n)Ho = {o, a^}\n)Ho, where a, a"^ are the standard bosonic annihilation 
and creation operators, the anti-commutation relation for harmonic oscillator gives its place to 
commutation relation for pseudoharmonic oscillator. 



5 Coherent states for pseudoharmonic oscillator: Barut- 
Girardello and Gilmore-Perelomov 

We are now able to produce two distinct classes of coherent states for the considered system. 
Usually, the Barut-Girardello type of coherent state is defined as the right eigenstate of the 
annihilation operator, i.e., 

M_\z)bg = z\z)bg- (44) 
To calculate the explicit form of \z)bg with the traditional method, we write 



oo 



\z)bg = ^Cn\n). (45) 



n=0 



Setting it in (14^ and making use of (1371) finally arrives us at the explicit form of coherent state 

T(s + 



^^)bg = Nbg{\z\')Y1 



n=0 



3n nl/2 
2/ 



_r(n + i)r(n + s + |)_ 



z"|n), (46) 



where the factor NbgH^I ) can be calculated by the normalization condition: {z\z) = 1, i.e. 

NBG{\zf)=(^oF,is + ^,\z\'?j (47) 

where oFi{a,x) denotes the hypergeometric function. 

On the other hand, referring the su{l, 1) symmetry algebra of the generators as in (H2l) . we 
are able to construct the Gilmore-Perelomov coherent state associated with the pseudoharmonic 
potential as follows [55] : 

\z)gp = I)(^)|0)=exp(eM+-rM_)|0) 

= exp(zM+)(l - |zn*^Oexp(-z*M_)|0), (48) 
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with z = ill tanh(^). Straightforward calculations leads one to the corresponding normalized 
Gilmore-Perelomov coherent state, i.e., 



r(n + s + f) 



,1/2 



r(n + i)r(s + f) 



\z)gp = Ngp{\z\'')J2 

n=0 

with the normalization factor 

Nap{\z'\)={l-\z\') 



z-\n), 



s_i3 

4 



(49) 



(50) 



Now, after finding out the two distinct classes of coherent states, we are ready to investigate 
some of the nonclassical properties which are of the most importance in the content of coherent 
states in quantum optics studies. But, before paying attention to this subject, we should inves- 
tigate the resolution of the identity for the above introduced states. If we became successful, 
the introduced states can be truly named as coherent states. 

Resolution of the identity: Any coherent state 1^) should satisfy the resolution of the identity 
which is defined as follows 



d^zw{\z\^)\z){z\ 



^|n)(n| 



(51) 



n=0 



where w(|2;p) is a positive weight function that we search for obtaining it. For the first class 
of states which is introduced in this section, substituting ( H^ in (1^ and using the change of 
variables z = re*^ and = x lead us to the following integral relation 



T{n + l)T{n + s + 

f(^+i) 



(52) 



where we have used the abbreviation wbg{x) = TrwBG{x)N'^Q{x). This is indeed an inverse 
moment problem which can be solved by the well-known methods such as Mellin transform 
technique |56j- However, using the definition of Meijer's G-function, it follows that [57j 



dx X 



k-l 



(3x 



ai. 



■ 1 l^m 1 l^r. 



■■■■> Oip 



(53) 



1 nr=ir(&, + fe)n?=ir(i-a,-fc) 
fi' n^=.+i r (1 - h, - k) r (a, + k) ■ 



Comparing equations (152|) and (153|) . it can be easily checked that the required weight function 
which can satisfy the integral equation (15T|) reads as: 



wbg{x) 



3 2^ 
2 ' 



7rr(s + 



-G 



4,0 
2,4 



X 



(0,3 + 1/2) 
(0,0,5 + 1/2,5 + 1/2) 



(54) 



where G is the well-known Meijer function and QFi{a;x) is the generalized hypergeometric 
function. 

About the Gilmore-Perelomov states, by substituting fH9|) in f l5T|) we will have: 



WGp{x)x^dx 



r(n 



r(n + s + |) 



(55) 
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where we have used wgp{x) = ttwgp{x)Nqp{x). Similarly, in this case it is easy to check that 
the required weight function has the following form 

where 2-^1(01, 02; a;) is the generalized hypergeometric function. 

A few words seem to be necessary about the limits of the integrals in (152|) and ( |55|) . Indeed, 
any quantum state belongs legally to a Hilbert space if its normalization constant satisfy the 
inequality < < 00. This led us to the unit disk and whole of the complex plane respectively 
for Barut-Girardello and Gilmore-Perelomov state. 



6 Nonclassical features of the introduced states 

In this section, we briefly discuss some of the criteria of nonclassicality which are usually used 
in the relevant literature [58], and will be helpful for investigating the nonclassicality exhibition 
of our introduced states in the continuation of the paper. To achieve this aim, we refer to 
the normal squeezing, amplitude-squared squeezing and sub-Poissonian statistics as the most 
popular nonclassicality signs. To mention the common feature of the above criteria we should 
indicate that if a state is squeezed or its quantum statistics is sub-Poissonian in a region of 
allowed space, the corresponding Glauber-Sudarshan P-function will not be positive definite, 
i.e., it has no classical analogue. But, we would like to imply that finding this function is usually 
a hard task to do. Altogether, each of the above effects is in fact a sufficient (not necessary) 
condition for a quantum state to belong to the family of nonclassical states. 

• Normal squeezing: In order to examine the quantum fluctuations of the quadratures of 
the field [59] we introduce the Hermitian operators 

M_+M+ M_-M+ 
X^ = , P^ = (57) 

The squeezing parameters can be defined as follows 

S-y= , ^ -1, 7 = ^1,^1, (58) 

where {A^/y = (7^) -(7)^. Using §7!) in f l58l) we obtain the following first-order squeezing 
parameters: 

Sx. = ^^{(Ml) + {M_M^) + {M^M_) + (Ml) 

- (M!)-(Mn-(M_)(M+)-(M+)(M_))-l (59) 



= ^^(-(M2) + (M.M+) + (M+M_)-(M^) 

+ (M!) + (M^)-(M_)(M+)-(M+)(M_))-1, (60) 
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where all expectation values can be straightforwardly obtained with respect to both classes 
of states. A state is squeezed in Xi or Pi if it satisfies the inequalities — 1 < Sxi < or 
— 1 < Sp-^ < 0, respectively. 

By comparing the two displayed diagrams in figure (1) it can be realized that the squeezing 
of Gilmore-Perelomov coherent states occurs only in Pi direction in all real space. Our 
calculations about the squeezing for Barut-Girardello coherent states in Xi and Pi, show 
that this parameter is approximately 10~^^ (= 0) for both operators. Thus, in view of 
squeezing property, it is deduced that the Barut-Girardello states associated with the 
considered system behave similar to the linear harmonic oscillator. 

Amplitude- squared squeezing: Following the proposal of Hillery [60], to evaluate this pa- 
rameter we define the Hermitian operators 

+Ml M'^ - Ml 
X2 = —— — P2 = — -. 61 

The squeezing conditions in X2 or P2 are respectively given by — 1 < 8x2 < or — 1 < 
Sp^ < 0, where 8x2 and Sp^ defined as follows: 

Se= , ^^^^' _-l, £ = X2,P2 (62) 

with the definition (A£)^ = — The explicit form of the latter parameters for 
second-order squeezing can also be obtained in a similar manner to what is done in 
equations (159!) and (!60l) for the normal squeezing. In figure 2 amplitude-squared squeezing 
has been plotted against real z for Gilmore-Perelomov coherent states, in X2 and P2 
direction. According to the diagrams have been placed in figure 2, it appears that, like 
the first-order squeezing, for this state again amplitude-squared squeezing occurs only in 
P2 direction. Also, this parameter for Barut-Girardello states in the two directions X2 
and P2 takes approximately the value of zero. This latter observation is also the same as 
the situation which is likely occurred in ordinary coherent states. 

Mandel parameter: To establish the quantum statistical properties and as another non- 
classicality sign, we calculate the Mandel parameter [61], which can be generalized as: 

Q = ^ + — ^± L - 1. 63 

The state for which (5 = 0, Q < and Q > respectively corresponds to the Poisso- 
nian (standard coherent states), sub-Poissonian (nonclassical states) and super- Poissonian 
(classical states) statistics. In figures 3 Mandel parameter has been plotted versus real z 
for Gilmore-Perelomov coherent states. As is shown, these states behave nonclassically 
features in a finite region of space. Our numerical calculations of this parameter for 
Barut-Girardello states takes the fixed value Q = —1. Therefore, their quantum statistics 
is such that Mandel parameter gets the lowest possible negative value, like the number 
states of harmonic oscillator which sometimes have been called as the most nonclassical 
states. 



11 



Adding the above observations, we may conclude that, from the point of view of squeezing 
properties the pseudoharmonic oscillator Barut-Girardello coherent state which is a particular 
superposition of number states (of pseudoharmonic oscillator) behaves like the canonical co- 
herent states \z) = exp(— |zp/2) Yl'^=o '^\^) ho- However, taking into account the quantum 
statistical properties of the field, the pseudoharmonic oscillator Barut-Girardello coherent state 
behaves like the number states of harmonic oscillator (note that \n)Ho 7^ 1^)) where \n) in the 
present paper is related to the eigenkets of pseudoharmonic oscillator potential). So, in this 
view, the latter state can also be named as the most nonclassical states, too. 



7 Summary and conclusion 

In summary, in this paper we obtained the eigenvalues and eigenfunctions of one-dimensional 
pseudoharmonic oscillator by the factorization method. The presented approach allows us to 
find out appropriate generalized creation and annihilation operators by which one can raise and 
lower arbitrary eigenkets of the considered system. Then, we have explored that the system 
under consideration possesses a symmetrical realization of a particular dynamical group; it is 
shown that the constructed ladder operators can constitute the generators of the su{l, 1) Lie 
algebra. Next, we have produced the explicit form of Barut-Girardello and Gilmore-Perelomov 
coherent states for the considered system. The resolution of the identity, as the most important 
task for establishing any state as a class of "coherent states" is estabhshed analytically and 
so the over-completeness relation for both states are demonstrated. Finally, we could high- 
lighted some of the nonclassical properties, such as first-order, amplitude-squared squeezing 
and sub-Poissonian statistics for the introduced states, numerically. Interestingly, we have 
also found that, in particular, considering squeezing properties, the pseudoharmonic oscillator 
Barut-Girardello coherent state behaves such as the canonical coherent states, however taking 
into account the quantum statistical properties, the latter state behaves like the number states 
of harmonic oscillator. Although, we have applied the method to pseudoharmonic oscillator in 
this paper, the approach can be used for other systems of interest, too. As another example, we 
have also done the same procedure to the well-known problem: a quantum particle in an infinite 
square well potential, successfully, with the new and novelty results, which will be appeared in 
near future. 
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Figure Captions: 

Fig. 1 Plot of squeezing parameters, Sx^ (up diagram) and Sp^ (down diagram) against 
real z with s = 1 for Gilmore-Perelomov coherent states. 

Fig. 2 Plot of amplitude-squared squeezing parameters, 8x2 (up diagram) and S'pj (down 
diagram) against real z with s = 1 for Gilmore-Perelomov coherent states. 

Fig. 3 Plot of Mandel parameter with s = 1 for Gilmore-Perelomov coherent states versus 
real z. 



15 



Sxi ,Sp^ 




z 



5 
4 
3 

Q 2 
1 

-1 



